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In  1954  Volume  51  of  the  "Notes  of  the  All-Union  Scientific 
Research  Institute  of  Hydraulic  Engineering"  contained  an  article 
which  examined  the  method  for  calculating  the  freezing  rate  of 
water-permeable  ground  by  a row  of  columns  after  frozen  ground 
cylinders  had  joined.  This  article  proposes  a method  for 
calculating  the  rate  at  which  water-permeable  ground  is  frozen 
by  a series  of  columns  before  the  frozen  ground  cylinders  join. 

1.  Statement  of  the  Problem 


The  problem  under  examination  can  be  formulated  as  follows. 


Let  there  be  given  a cofferdam  of  length  L and  thickness  H (Figure 
1).  The  cofferdam  is  built  on  a base  which  is  permeable  to  water  and 
heat.  Water  filters  through  the  cofferdam  under  the  action  of  a 
constant  difference  in  head  between  the  lateral  boundaries  of  the 
cofferdam.  Water  does  not  filter  around  the  cofferdam.  The  freezing 
columns  with  radius  r^,  height  h and  total  number  L are  arranged  on 

the  cofferdam  on  a plane  which  is  parallel  to  the  lateral  (head) 
surfaces  at  uniform  distances  S from  each  other.  Before  the  columns 
are  implanted  in  the  ground,  the  filtration  flow  under  natural 
(normal)  conditions  has  speed  v^  and  temperature  T.  The  ground 

freezes  at  temperature  t^.  The  thermophysical  parameters  of  the  frozen 


and  unfrozen  ground  are  known.  It  is  necessary  to  find  a dependency 
which  determines  the  ground  freezing  rate  before  the  cylinders  join. 
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Figure  1.  a.  Plan  of  Cofferdam;  b,  Transverse  Cross-Section 
of  Cofferdam. 
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2.  Heat  Balance  Equation 


The  heat  balance  equation  of  the  frozen  ground  looks  as  follows: 


^f 


! n 


where  Q^.  is  the  amount  of  heat  evolved  per  unit  of  time  when  water  freezes 
on  the  surface  of  the  frozen  ground  cylinders; 


CESl  AVAIl/i^it 


is  the  amount  of  heat  which  is  given  off  per  unit  of  time  by  the 
freezing  columns  (the  thermal  absorption  by  the  freezing  columns); 

is  the  amount  of  heat  which  the  filtration  flow  gives  off 
per  unit  of  time  to  the  frozen  ground  cylinders; 

Qg  is  the  amount  of  heat  which  the  frozen  ground  loses  per  unit  of 
time  (cooling  of  the  frozen  ground  cylinders). 


The  thermal  flow  Q^,  as  is  demonstrated  in  examining  the  freezing 

of  water-permeable  ground  by  a single  column  [1]  can  be  taken  into 
account  by  increasing  the  calculated  value  of  the  latent  ice  formation 
heat  by  the  value 


Ur 


g' 


(2) 


where  is  the  thermal  conductivity  factor  of  the  frozen  ground; 
c^  is  the  thermal  capacity  of  the  frozen  ground; 

is  the  specific  weight  of  the  frozen  ground; 
h is  the  height  of  the  frozen  ground  cylinder. 

Therefore  thermal  balance  equation  (1)  can  be  represented  in  the 
following  simpler  form: 


Qf  = Qc  - Q,i-  ( 

The  thermal  flux  from  the  latent  ice  formation  heat  is  equal  to: 


i-n 


Here  R is  the  radius  of  the  frozen  ground  cylinder; 

j is  the  latent  ice  formation  heat  of  a unit  of  ground  volume; 


r is  time. 


Allowing  for  the  thermal 

write 


flux 


instead  of  (4)  it  is  necessary  to 


where 


(01 
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The  absorption  of  heat  by  the  freezing  columns  changes  as  the 
ground  freezes.  The  law  governing  the  change  in  heat  absorption  is 
determined  by  the  characteristic  of  the  refrigeration  device  and  by 
the  thermal  resistance  of  the  frozen  ground. 


.3 


The  means  for  determining  this  dependency  are  analyzed  in 
the  technical  literature  [2]:  here  we  will  write  this  dependency  with 
regard  to  the  frozen  ground  cylinders  in  the  following  general  form: 


Qc  = f(R).  7 

This  expression  may  be  regarded  as  a known  dependency;  in  doing  this  1 

we  keep  in  mind  the  fact  that  by  properly  selecting  the  refrigeration 
devices  this  dependency  can  be  given  any  desired  form. 

) 

If  the  freezing  is  conducted  at  a constant  brine  temperature  9,  ^ 

then 


'•u 

It  remains  for  us  to  determine  the  value  of  the  thermal  influx 
from  the  filtration  flow  The  determination  of  this  component  of 

heat  balance  equation  (3)  is  the  basic  difficulty  in  solving  the  stated 
problem. 

3.  Calculating  the  Thermal  Influx  from  the  Unfrozen  Water-Permeable 
Ground  to  the  Frozen  Ground  Cylinders 

We  will  determine  the  value  of  the  thermal  flux  from  the  unfrozen 
ground  to  the  frozen  ground  cylinders  on  the  basis  of  the  following 
equation: 

where  c^  is  the  thermal  capacity  of  the  water; 
is  the  specific  weight  of  the  water; 

I is  the  number  of  columns; 

Vg  is  the  speed  of  the  filtration  flow  in  the  cofferdam  in  the 
area  where  the  frozen  ground  cylinders  do  not  cover  the  entire  cross-section; 

S is  the  distance  between  the  axes  of  the  freezing  columns; 

T is  the  temperature  of  the  filtration  flow  at  the  approach  to  the 
frozen  ground  cylinders  (equal  to  the  temperature  of  the  filtration  flow 
under  natural  conditions); 

t^^^  is  the  average  temperature  of  the  filtration  flow  after  passing 
through  the  frozen  ground  screen. 
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In  order  to  make  subsequent  analysis  convenient,  we  will  introduce 
two  parameters; 

_ V. 

' ^ I 

"•b" 


(8) 


where  is  the  speed  of  the  filtration  flow  before  freezing  begins  and 


\V~J1 ^ut  . 
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(9) 


where  as  before  t^  is  the  freezing  temperature  of  the  ground. 


Then  the  equation  for  the  thermal  flux  acquires  the  following 


form: 


where 


:J.-  c V S/.'I-j  { F-t,,). 

w w b 


(10) 


(11) 


In  order  to  use  equation  (10)  it  is  necessary  to  know  the  values 
of  V and  At. 


In  the  problem  which  we  are  examining  (a  row  of  freezing  columns, 
without  filtration  around  the  cofferdam)  as  the  frozen  ground  cylinders 
increase  in  size,  the  rate  of  flow  through  the  cofferdam  decreases, 
and  therefore  the  speed  of  the  filtration  flow  v^  and  the  value  of  v 

decrease.  This  problem  was  studied  by  G.  S.  Shadrin  who  used  the 
graphic  method  of  constructing  motion  grids  [3].  An  analytic  examination 
of  this  problem  has  led  us  to  the  following  relationship  [4]  which  determines 
the  desired  parameter: 


where 


■ 


(12) 


(13) 


(14) 


( 15) 


A graph  of  the  dependency  v = f(r,  is  presented  in  Figure  2. 
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Figure  2.  Graph  for  Determining  the  Relative  Change  in 
Filtration  Flow  Through  the  Cofferdam  When  the  Frozen  Ground 
Cylinders  Increase  in  Size. 

Let  us  know  determine  At.  For  this  purpose  we  will  use  the  solution 
given  by  heat  transmission  theory  for  the  problem  of  wall  cooling 
[5].  The  problem  is  formulated  as  follows. 

An  infinite,  homogeneous  wall  of  thickness  2X  at  the  initial 
time  (t  = 0)  has  the  same  temperature  T at  all  points,  and  then  (t  > 0) 
both  sides  of  its  surface  are  e.xposed  to  a temperature  of  t^.  The 

temperature  conductivity  factor  of  the  wall  a is  known. 

The  solution  to  this  problem  has  the  following  form: 

(16) 

where 


A graph  of  dependency  (16)  is  presented  in  Figure  3,  the  lower  left 
quadrant  (see  glued-in  piece  at  the  end  of  the  book). 

We  will  show  how  to  use  this  solution  in  our  case. 

We  will  write  the  differential  equation  of  heat  conductivity  in 
water-permeable  ground  (the  arrangement  of  the  coordinate  axes  is  given 
in  Figure  4) ; 

:-r  ■:  , 

(IS) 
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Ignoring  the  value  of  the  thermal  flux  which  arises  due  to  physical 
heat  conductivity  and  is  directed  along  axis  Y and  the  component  of 
induced  convection  which  is  directed  along  X,  we  will  convert  equation 
(18)  into  the  following  form: 


■jy 


where 


a 


(19) 


(20) 


is  the  thermal  conductivity  factor  of  unfrozen  ground. 


Figure 

Equation  (19)  can  be  regarded  as  a one-dimensional  thermal  conductivity 
equation  in  which  the  time  variable  is  replaced  by  the  y coordinate 
variable. 

Thus,  the  task  of  cooling  a filtration  flow  which  occurs  iu  a gap 
2X  thick  amounts  to  the  task  of  cooling  a wall,  the  solution  to  which 
is  known  (equation  (16)).  However,  in  the  case  which  is  of  interest 
to  us,  cooling  a filtration  flow  when  it  passes  through  a frozen  ground 
stream,  we  are  dealing  with  a gap  of  non-constant  thickness. 

We  will  show  that  the  above- indicated  solution  can  also  be  applied 
to  our  problem. 
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i 

t 

1 


^ n I*' 


As  is  known,  the  temperature  distribution  in  an  infinite  planar- 
-parallel  wall  is  determined  by  integrating  the  differential  equation 


d-  ox  ' (21) 

If  we  place  the  origin  of  the  coordinate  axis  in  the  middle  of  the 
wall  (x  = 0)  and  designate  the  wall  thickness  as  2X,  then  under 
boundary  conditions 


n.v,  0,^7^  (22) 

'•—.V,  T)  .ri  -..V.  T)=.r„;  (23) 

the  solution  to  equation  (21),  as  is  well  known,  will  be  as  follows  [5]: 


where 


(24) 


(25) 


(26) 


"'I 


(27) 


(28) 


For  any  moment  in  time  (t  = t^)  equation  (24)  can  be  regarded  as 

an  initial  condition  (instead  of  (22))  with  equation  (21)  integrated. 

This  initial  condition  is  characterized  completely  by  the  value  of 
criterion  Fo;  in  this  case,  from  the  viewpoint  of  the  subsequent 
process  it  is  completely  irrelevant  what  the  conditions  are  under  which 
the  process  of  heat  transmission  has  occurred  during  the  preceding 
period.  The  course  of  the  wall's  temperature,  as  before,  will  be 
described  by  equation  (24),  and  the  speed  of  increase  of  Fo  will  depend 
on  the  wall  thickness  at  the  given  time.  We  will  consider  that  the 
change  in  wall  thickness  and,  consequently,  the  change  in  the  area  of 
the  filtration  flow's  cross-section  per  se  does  not  cause  any  temperature 
changes  at  compatible  points,  i.e.,  we  will  assume  that  when  the  stream 
of  the  filtration  flow  is  instantaneously  compressed  or  expanded, 
dependency  t = f('r^)  will  remain  unchanged. 

In  this  case 

1 29) 
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And,  consequently, 


(.30) 


Therefore  the  value  of  criterion  Fo  which  has  to  be  substituted  into 
the  solution  of  (24)  is  determined  by  the  integral: 


I'lj  - n 


(51) 


In  order  to  obtain  this  integral,  we  will  express  dt  through  X, 
and  in  doing  so  we  note  that 


From  Figure  4 it  follows 


(32) 


y^  R-  I R--- 


by  differentiation,  we  find: 


(35) 


Keeping  in  mind  the  fact  that  the  average  velocity  of  the  filtration 
flow  is  inversely  proportional  to  the  area  of  the  cross-section,  we  can 
write : 


By  substituting  (32)  into  (31)  and  taking  |33)  and  (34)  into 
account,  we  find: 


Bv  integrating  within  the  limits  of  X = t to  ^ doubling  t!ic 

result  (since  the  limits  cover  only  half  of  the  length  of  the  slit  I and 
bv  substituting  expression  1 13),  we  find  the  value  ot  the  Fourier 
criterion  after  the  filtration  flow  passes  through  the  frozen  ground 
screen : 


i 


I 

i 


I 


T » • 

uL. 

and,  keeping  (12)  in  mind,  we  finally  obtain: 

(■/,  -I, 


where 


(35) 


(36) 


R'- 


(37) 


--  see  equation  (15). 

The  results  of  calculating  the  Fourier  criterion  in  terms  of 
dependency  (35)  are  also  presented  in  Figure  3.  By  knowing  x.  e !ind 
K,  this  dependency  makes  it  possible  to  find  At. 

Thus,  by  using  the  graph  in  Figure  2 it  is  possible  to  determine 
the  values  of  parameter  v,  and  from  the  nomogram  in  Figure  3 it  is 
possible  to  find  the  values  of  parameter  At.  This  is  sufficient  to 
calculate  the  magnitude  of  the  heat  influx  from  the  unfrozen  water- 
-permeable  ground  according  to  formula  (10).  Below  is  given  a 
corresponding  sample  calculation  for  the  purpose  of  illustration. 

4.  Sample  Calculation  of  Heat  Influx  From  Unfrozen  Water-Permeable 
Ground 

We  will  determine  the  thermal  influx  to  one  of  the  frozen  ground 
cylinders  which  is  in  a row  with  others. 

The  following  are  given:  V.  = 1 — — ; 

* * 2 m-hr-deg 


, kcal 

c = 1 — j — ; 
w m-deg 


m 


V,,  = 0.04  m/hr;  T = ~°C;  t,^  = 0°C;  S = 1 m;  H = 3 m;  I = 1;  h = 1 m. 
b 0 

From  equations  (20),  (14),  (37),  and  (11),  we  find: 


7 

H_ 

S 


I I 1 • 

W'W 


• hr 


kcal  hr. 
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According  to  computation  equation  (10) : ” - ^ 

\o' 

The  rest  of  the  computations  are  summarized  in  the  table.  For  the 
sake  of  comparison  the  last  graph  in  the  table  shows  the  results  of 
calculating  the  heat  influx  to  a single  frozen  ground  cylinder,  with 
these  computations  carried  out  according  to  formula  of  B.  V. 
Proskuryakov  [6] : 


k'T-r  .)  1 . 


(38) 


TABLE  1. 


I 

I 


Commas  indicate  decimal  points. 

From  Table  1 it  is  evident  that  the  magnitude  of  the  heat  influx 
from  the  unfrozen  ground  to  a cylinder  which  is  in  a row  with  another 
reaches  its  peak  value  when  e ^ 0.8,  and  then  rapidly  decreases. 

A comparison  with  the  data  for  a single  cylinder  shows  that  in  the  cited 
case  conditions  for  ground  freezing  by  a row  of  columns  are  more  favorable. 
The  corresponding  calculations,  however,  show  convincingly  that  under 
these  circumstances,  for  instance  at  higher  filtration  flow  speeds  or 
when  the  cofferdam  is  quite  thick,  cases  are  possible  where  the  thermal 
influx  to  the  frozen  ground  cylinders  in  a row  with  other  cylinders 

is  greater  than  the  flow  to  a single  cylinder. 

5.  Solution  to  the  Heat  Balance  Equation 


We  will  substitute  into  heat  balance  equation  (3)  the  expressions 
of  its  components  (5)  and  (10),  and  by  taking  (13)  into  account  we 
find; 


••'c 


\ z- 


The  solution  to  the  posed  problem,  determining  the  ground  freezing 
time  before  the  frozen  ground  cylinders  ioin,  is  found  by  integr.it  ing 
this  latter  equation  within  the  limits  from  = dr^^.'S  to  = 1 and  from 

• =0,  to  - = • 
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(39) 


• I 
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The  sub-integral  function  is  complex  in  form,  and  therefore  the 
value  of  the  integral  must  be  calculated  by  the  graphic  or  tabular 
method,  and  then  determining  the  freezing  time  according  to  formula 
(39)  presents  no  problem. 

6,  Sample  of  Calculating  Freezing  Rate 


We  will  determine  the  period  of  ground  freezing  which  occurs  before 
the  frozen  ground  cylinder  join  under  the  following  initial  conditions: 


length  of  cofferdam 

thickness  of  cofferdam 

height  of  column 

distance  between  column  axes  

radius  of  freezing  columns 

speed  of  filtration  flow 

temperature  of  filtration  flow  

latent  ice  formation  heat  per  unit  of 

ground  volume 

thermal  capacity  of  water  

specific  weight  of  water 

thermal  conductivity  factor  of 
unfrozen  ground 

thermal  conductivity  factor  of 
frozen  ground 

thermal  capacity  of  frozen  ground.... 

specific  weight  of  frozen  ground 

We  find 


L = 100  m 
H = 6 m 
h = 8 m 
S = 1 m 
r^  = 0.05  m 

V,  = 0.08  m/hr 
D 

T = 6“C 

o = 24.000  kcal/m'^ 
c^  = 1 kcal/kg/hr 

Y = 1000  kg/m'^ 

X.,  = 1 kcal/m-hr-deg 

\^  = 2 kcal/m’ hr* deg 
c^  = 0.34  kcal/kg°C 
= 16.00  kg/m"^. 


1 ■ lOXO 


M‘  hr 


w'w 

4i;  4 • O.COl 


0,08  . 1 

' S 1 


--0.0'); 


, L '.Oit  . 

\ -c  S/!/i'i  ! r.l-  1 lOOO- 1 • 100  ■ II.QX  h kcal  hr. 
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The  refrigeration  device  ensures  the  constant  brine  temperature 
9 = -20°C;  according  to  formula  (7)  and  taking  into  account  the  fact 
2^ 

that  R = ^ , we  will  determine  how  much  heat  the  freezing  columns 
absorb : 


'x'c 


. 100  ■2-20 


In 


2 ■ 0,05 ) 


201  • in- 

In  ( 10s ) 


The  calculated  value  of  the  latent  ice  formation  heat  is  determined 
from  equations  (2)  and  (6) : 


<?.  2'!000-- 
•1  ; In 


C,.j  t • l60 

4.r  . 


■0  2 4<>'j0-f-27 


111 


Qc 


The  rest  of  the  calculations  are  summarized  in  Table  2,  and  in  this 
process,  in  accordance  with  equation  (39)  we  utilize  the  following 
expression: 


.V-A:  2 ‘ 


-I26U:' 


table  2, 


1 . 
1 ' 

-i  . 

kcal 

1 = 

Found 

Trom 

■fj-.V 

kcal 

'c  fi 
kcal 

'e  -Ti 

hr 

j 

' kcal 

' ■ M 

Accord- 
ing to  , 

fj^rmiiia  nrs 

1 

1 

hr 

figure 

r 

-TT 

”KT 

nr 

kcal 

iUiJJiUid 

(39),  hn 

K- 

".1 

2!,010' 

0.:O 

i o.ir- 

!3700 

:4s.ii!) 

; U.'l  ill  ' 

: r.'i..iii 

i'.i.i 

KSMO". 

! 0.100 

,'OSOo 

17 -'.uo 

J,.5  1(1' 

' AS'  li  , '.1 

io.‘, 

■I'.l 

14.'  ’.0" 

0,97 

O.i'lj 

7.ii!i)0 

7:ai  'I 

1 

'7l<  *m 

• ■ . 1 

’ .-n 

[ar, 

O.t 

125- 10- 

o,:'{ 

1 

■'0500 

3'5o; 

' 1 J 1 Cl  ■ ’ 

1 LT370 

, , , 

ja« 

(t.l 

I ’->•  :o- 

C'l 

o.'joo 

lonoo 

107C0 

j5G  111 

i ■J7i;i'0 

JOii.sy 

! 

110- 10- 

0,'4. 

0 090 

•111'  00 

'li'iO 

1 i:-  ; ■ in 

1 'j>i..:ii 

' 1 

.UVb'-. 

i'Vfw. 

107  I o- 

OS' 

. 0 .41 

111.50(1 

1 ■_  j'-nO 

' 

Commas  indicate  decimal  points. 

Thus,  under  the  given  conditions  joining  will  not  occur; 
the  frozen  ground  cylinders  will  cease  to  grow  as  early  as  e — O.bd. 

In  the  given  case  it  is  advisable  to  shift  to  freezing  at  a lower 
brine  temperature.  Below  Table  3 shows  a calculation  of  freezing  time 
which  demonstrates  that  if  we  assume  that  the  brine  temperature  is 
9 = -40“C,  then  joining  will  occur,  for  instance,  on  the  loth  day. 


-13- 


i 


BEST 


# f I 

^ A i 


hu:  ^ ^ i 


TABLE  3. 


1 , 

( 

i Found 
1 from 
jF^gure, 

- 

if 

Found 

from 

^gure 

ij  • .\  - a7. 

kcal 

~Tir 

tj  - rj  , 

kcal 

■ 

0 -'5 

hr 

kcal 

. a* 

kcal  According 
— 3"  CD  formula 
' (391.  hrs 

hrs 

o.n.'i 

457CO 

536300 

0,37 

10-' 

.39700  ! 

1.85 

1.5.3 

J i).44 

II.I6U 

rjhoii 

.30.5.00 

0.9S 

;o"' 

' 33900 

).:o  } 

6,0.8 

I'.'.A 

73000 

J170U0 

1 <4 

10  ■ 

: .3 If. 40 

- , . 1 
. . ih  * 

13.41 

■ 

O.Jiiy 

.SfiSOO 

163500 

!.05 

•i)-' 

307c*  > 

118  i 

25.21 

nvl 

1)  .'90 

101300 

122700 

4.8 

lo-‘ 

' 300.30 

18.12  ! 

43,.i3 

■ 

y.3o'. 

i 1 6000 

yoooo 

7.73 

10  •' 

i ■;'95.30 

25,95  1 

72.28 

■ 0.7  7 

0,  i'5o 

1.36000 

58000 

13.8 

10  ■ 

292.;o 

'064  ' 

1.2,92 

D.tiZ 

0.690 

Ki^OOO 

20000 

45 

lO"' 

i 28970 1 

164.0  1 

2£6.<J2 

’ '1,47 

‘L'^yo 

1 (X16IX) 

17-100 

2.3,8 

10-' 

I 28.800  1 

S6..2  * 

373.22 

■ 0,  ,.{ 

0.095 

I'JtiMOO 

50000 

5.88 

10-' 

■ 2.S7.'.0  • 
, 1 

2 1. 3.3 

0 : '7 

Commas  indicate  decimal  points, 
loining  Conditions 

r the  cylinders  to  join  into  a single  frozen 
of  the  fact  that  the  heat  absorbed  by  the  columns 
; must  exceed  the  thermal  influx  from  the  filtration 

iwing  condition  must  be  fulfilled: 


) 


■J:  fT  ' ■ O / 

Assume  that  it  is  necessary  to  determine  the  distance  between  the 
freezing  columns  at  which  joining  of  the  frozen  ground  cylinders 
will  be  ensured,  if  it  is  known  that: 

length  of  cofferdam L = 50  m 

thickness  of  cofferdam H = 8 m 

height  of  freezing  columns h = 10  m 

radius  of  freezing  columns r^  = 0.05  m 

speed  of  filtration  flow = 0.05  m/hr 

temperature  of  filtration  flow T = 5°C 

thermal  capacity  of  water ~ ^ kcal/kg-deg 

specific  weight  of  water = 1 kg/m"^ 

thermal  conductivity  factor  of 

frozen  ground = 2 kcal/ra* hr* deg 

thermal  conductivity  factor  of 

unfrozen  ground = 1 kcal/m- hr' deg 

brine  temperature  in  columns .- . 6 = -22°C 

freezing  temperature  of  ground t^  = 0°C. 

According  to  (11)  and  bearing  in  mind  the  fact  that  I = L/S, 
we  can  write 

• iO  ■ 50  ■ 0/>5  ■ ><cal/hi. 

If  we  assume  that  S = 1.5  m,  then  according  to  e.xpression  (57)  and 
(14)  [with  allowance  for  (20)]: 


'X  • 1.0 


1 • 1UOO.U,05.  1.^ 


n „ 

From  Figure  5 we  determine  that  the  maximum  heat  influx  of  the  filtration 
flow  occurs  when  e — 0.94. 

From  Figures  2 and  3 we  find  that  when  z s 0.94: 

V = 0.48  and  At  = 0.86,  and  consequently  according  to  (10): 

' fi  0,.r-  .!  I, S' • .ii.i,'  kcal/hr. 

On  the  other  hand,  the  heat  absorbed  by  the  columns  is  equal  to: 


• "0  ■ 2.0  • . 


kcal/  hr. 
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Thus,  it  has  been  found  that  when  S = 1.5  m ^ consequently 

joining  does  not  occur;  therefore  it  is  necessary  to  set  the 
columns  closer  together.  If  we  assume  that  S = 0.8  m and  repeat  the 
course  of  the  above-cited  calculation,  we  find  that  the  position 

If  <a  1 

c = 0.92  will  be  most  difficult  and  in  this  case:  Q...  = 75,000  — ^ 

hr 

and  = 86,500  Therefore  this  version  is  technically  feasible 

since  condition  (40)  is  maintained. 


Figure  5.  Graph  for  Determining  tlie  Value  of  Parameter 
e at  Which  the  Thermal  Influx  From  the  Unfrozen  Water- 
-Permeable  Ground  to  the  Frozen  Ground  Columns  Reaches 
Its  Peak. 

Cone  I us i ons 

In  this  article  we  have  examined  the  problems  of  making  a 
thermal  engineering  calculation  of  water-permeable  ground  being 
frozen  by  a row  of  columns  before  the  frozen  ground  cylinders 
join  where  the ^filtration  flow  arises  through  the  cofferdam 
under  the  action  of  acconstant  head. 

As  a consequence  a method  corresponding  to  this  case  has  been 
developed  for  calculating  the  period  of  freezing.  In  this  process 
it  was  found  that  as  the  frozen  ground  cylinders  increased  in  size, 
the  thermal  influx  from  the  unfrozen  ground  first  increased  and  then, 
due  to  the  reduction  in  the  flow  rate  of  the  filtration  flow,  decreased. 
Therefore  the  most  difficult  period  of  freezing  occurs  not  directly 
prior  to  the  joining  of  the  frozen  ground  cylinders,  but 
somewhat  beforehand.  If  by  the  time  the  maximum  thermal  influx  from 


-l(w 


the  filtration  flow  occurs  the  frozen  ground  cylinders  are  still 
growing  in  size,  then  it  is  considered  that  their  joining 
into  a unified  frozen  ground  mass  is  ensured. 

It  should  be  noted  that  the  method  proposed  in  this  work  for  a 
thermal  calculation  of  the  cooling  of  a filtration  flow  flowing 
between  two  frozen  ground  cylinders  may  also  be  used  to  solve  a 
number  of  other  thermal  problems  in  which  the  temperature  conductivity 
of  a body  and  its  dimensions  vary  according  to  a known  law  (for 
instance,  in  the  problems  of  water  course  thermics).  In  this 
process,  however,  it  must  be  kept  in  mind  that  this  method  is 
correct  only  when  there  are  constant  temperatures  at  the  boundaries 
and  if  the  physics  of  the  event  under  consideration  makes  it  possible 
to  consider  that  the  change  in  wall  thickness  (or  water  course  depth, 
diameter  of  the  penstock,  etc.)  does  not  per  se  cause  temperature 
changes  at  compatible  points. 
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